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Abstract 

A simple quark-diquark model of baryons with direct and exchange interactions has been con- 
structed. Spectrum and form factors have been calculated and compared with experimental data. 
Advantages and disadvantages of the model are discussed. 
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The notion of diquark is as old as the quark model itself. Gell-Mann [l| mentioned the 
)ossibility of diquarks in his original paper on quarks. Soon afterwards, Ida and Kobayashi 
3] and Lichtenberg and Tassie p introduced effective degrees of freedom of diquarks in order 
to describe baryons as composed of a constituent diquark and quark. Since its introduction, 
many articles have been written on this subject 4 1 up to the most recent ones P|. 

Different phenomenological indications for diquark correlations have been collected dur- 
ing the years, such as some regularities in hadron spectroscopy, the AI =_| rule in weak 
non-leptonic decays 3], some r^ularities in parton distribution functions [3] and in spin- 
dependent structure functions [2|- Finally, although the phenomenon of color superconduc- 
tivity 0] in quark dense matter can not be considered an argument in support of diquarks 
in the vacuum, it is nevertheless of interest since it stresses the important role of Cooper 
pairs of color superconductivity, which are color antitriplet, flavour antisymmetric, scalar 
diquarks. 

The introduction of diquarks in hadronic physics has some similarities to that of correlated 
pairs in condensed matter physics (superconductivity Q]) and in nuclear phvsics (interacting 
Doson model jl^l), where effective bosons emerge from pairs of electrons ^| and nucleons 



12| respectively. However, while the origin of correlated electron and nucleon pairs is clear 
(the electron-phonon interaction in condensed matter physics and the short-range pairing 
interaction in nuclear physics), the microscopic origin of diquarks is not completely clear 
and its connection with the fundamental theory (QCD) not fully understood, apart from 
the cold asymptotic high dense baryon case, in which the quarks form a Fermi surface and 
perturbative gluon interactions support the existence of a diquark BCS state known as color 
superconductor isj. We can only say that in perturbative QCD the color antitriplet, flavor 
a„t.yn,n,etnc =ca,a. channel i. faw.ed by one glnon exchange Q and in non-pertn.bat.ve 
regime by instanton interactions [1^. Regarding the hadron spectrum we are interested in 
the role of diquark correlations in non perturbative QCD. In this respect it is interesting 
the discussion of non-perturbative short-range, spin and flavour dependent correlations in 
hadrons, given in Ref.^^ considering an Instanton Liquid Model and the comparison of 
some of these results with LQCD ones. However, the question remains open and other 
lattice QCD calculations are needed in order to understand it fully. 

Nonetheless, as in nuclear physics, one may attempt to correlate the data in terms of a 
phenomenological model. In this article, we address this question by formulating a quark- 
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diquark model with explicit interactions, in particular with a direct and an exchange in- 
teraction. We then study the spectrum which emerges from this model and we start to 
calculate the form factors which have been measured, or will be measured at dedicated 
facilities (TJNAF,MAMI, ...)• 

We think of a diquark as two correlated quarks treated as a point-like object, thought 
this is a rough approximation of an extended effective boson degree of freedom. 

We assume that baryons are composed of a constituent quark, q, and a constituent 
diquark, = Q^. We consider only light baryons, composed of (u, d, s) quarks, with 
internal group SUs{2) SUf{3) SUc{3). Using the conventional notation of denoting spin 
by its value and flavor and color by the dimension of the representation, the quark has spin 
s = |, F = 3 and C = 3. The diquark must be 3 of SUc{S) since the total hadron must be 
colorless. This limits the possible SUsf{Q) representations for the diquark |^ to be only the 
21 of SUsf{6), which is symmetric and contains = 0, F12 = 3 and S12 = 1, -F12 = 6 . This 
is because we think of the diquark as two correlated quarks in an antisymmetric non-excited 
state. 

As one can see from the simple multiplication of the Young diagrams associated with two 
fundamental representations of SUsf{6) 

□ ® □ = m © (1) 

6 ® 6 = 21 © 15 , (2) 

by keeping only the representation 21 of SUsf{Q), as in the diquark case, one is deleting in 
baryons those states obtained by combining the representation 15 with that of the remaining 
quark 6, i.e. a 70 and 20. 

If one treats only non-strange baryons, the quark is in the representation 4 of the Wigner 
SUst{4:), with s = I and t = ^, and the diquark is the representations with spin S12 = 0, 
isospin ti2 = 0, and spin su = 1, isospin tu = 1, i.e. the symmetric representation 10 of 
SUst{4:) D SUs{2) SUt{2). The situation for the internal degrees of freedom is summarized 
in Table I. 

The two diquark configurations S12 = . F i? = 3 and S12 = 1, F12 = 6 are split, by 
(among other things) color magnetic forces as shown in Fig. 1. It means that we have 
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TABLE I: Combinations of possible irreducible representations of S'f7c(3), SUs{2) and SUf{3) for 
the quark, q, and the diquark, Q"^. The right-hand side of the Table also shows the combinations 
of the possible SUs{2) and SUt{2) representations if only non-strange baryons are involved. 

SUc{3) SUsf{6) SUs{2) SUf{3) SUst{4) SUs{2) SUt{2) 
q 3 6 s = i 3 4 s = ii=i 

Q2 3 21 si2 = 3 10 si2 = ti2 = 

Sl2 = 1 6 Sl2 = 1 tl2 = 1 



two possible diquark configurations: a scalar diquark (si2 = 0, F12 — 3, C — 3) and a 
vector one (si2 = 1, -F12 — 6, C — 3). The scalar diquark (the "good diquark" in Wilczek 
and Jaffe's terminology) is favoured since it is at lower energy, and will be the dominant 
configuration in the more stable states. 



_S12 = 1, ^^12 = 6 



C = 3 



A 



sio=0, F12 = 3 



FIG. 1: Schematic picture of the split, between the two diquark configurations S12 = 0, Fu = 3 
and S12 = 1, -F12 = 6 . 



We call the splitting between the two diquark configurations A and parametrize it as 

A = (S + Cdo) , (3) 

i.e. as a constant B which acts equally in all states with S12 = 1 plus a contact interaction 
which acts only on the ground state and of strength C. The contact interaction emphasizes 
the different role that the ground state has in a quark-diquark picture as compared with all 
other states. 

Here, we consider a quark-diquark configuration in which the two constituents are sepa- 
rated by a distance r. We use a potential picture and we introduce a direct and an exchange 
quark-diquark interaction. For the direct term, we consider a Coulomb-like plus a linear 
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confining interaction 

Vd^r{r) = + f3r . (4) 

The importance of tlie Coulomb-like interaction was emphasized long ago by Lipkin 
A simple mechanism that generates a Coulomb-like interaction is one-gluon exchange. A 
natural candidate for the confinement term is a linear one, as obtained in lattice QCD 
calculations and other considerations jisl . 

An exchange interaction is also needed, as emphasized by Lichtenberg This is indeed 
the crucial ingredient of a quark-diquark description of baryons. We consider 

Ve. = {-iy+'2Ae-'''[su ■ + h2 ■ h + 25^2 ■ sl,h2 ■ h]. (5) 

Here Si2, S3 and ti2 and ^3 are the spin and the isospin operators of the di quark and the quark 
respectively. The sign (—1)' reflects the angular momentum dependence of the exchange 
interaction. 

A feature of the present approach, as compared with previous ones, is that we attempt 
a simultaneous description of all states; indeed, in previous approaches the ground state 
has usually been excluded. It is certainly the case that diquark correlations are particularly 
important for high ^-states which have a more pronounced string-like behaviour. How- 
ever, if one is interested in the overall spectroscopy and also in the form factors (elastic and 
inelastic) that have been measured or will be measured (TJLAB, MAMI, ...) the ground 
state must be included in this description. 

Another important observation is that, in the quark-diquark model, SUsf{Q) is badly 
broken and thus a classification of the spin-flavor wave functions in terms of SUsf{Q) is 
inappropriate. The appropriate classification is that in terms of SUs{2) (g) SUf{3) 

I S12, F12; ^, 3; S, F) , (6) 

where S12 and F12 are, respectively, the spin and the flavor of the diquark, |, 3 are the spin 
and flavor of the quark and S, F are the total spin and flavor obtained by coupling the spin 
and flavor of the diquark with those of the quark. 

If only non-strange baryons are considered, one can use the isospin of the diquark, tu, 
and quark, 1/2, and the total isospin, T. The possible states are: 
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Sl2 = 0, ti2 
I Si2 = 1, ti2 
I = 1, tl2 
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(7) 
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The total wave functions are combinations of the spin-flavor wave functions with the 
radial and orbital wave functions. These are very simple and straigthforward since the 
main advantage of the quark-diquark model is to reduce the baryon problem (a three-body 
problem) to a two-body problem. The spatial part of the wave functions is, for the central 
interaction of Eq. (4), 

^n,«,™(rO = Rn,iir) YiA^,^) , (8) 

where the radial wave function Rn^i{r) can be obtained by solving the radial equation. 
Although the numerical solution poses no problem, we prefer to exploit here the special 
nature of the interaction. For a purely Coulomb-like interaction the problem is analytically 
solvable. The solution is trivial, with eigenvalues 



n 



1, 2 



(9) 



Here m is the reduced mass of the diquark-quark configuration and n the principal quantum 
number. The eigenfunctions are the usual Coulomb functions 



Rr,Jr) 



(^-^-^)'(2^)^-(2,r)^e-e-.(2.r). 



(10) 



2n[{n + l)\]^ 

n 

where for the associated Laguerre polynomials we have used the notation of Ref. [21J] and 
9 = —■ 

We treat all other interactions as perturbations. We begin with the linear term. The 
matrix elements oi (3 r can be evaluated analytically, with the result 



Pr[Rn,i{r)fr''dr 



--^-^[3.^-Z(/ + l)]. 



(11) 
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This perturbative estimate is only valid for small n and /. (For large n and /, the radial 
equation must be solved numerically.) Combining (9) with (11), we can write the energy 
eigenvalues as 

Next comes the exchange interaction of Eq. (5). The spin-isospin part is obviously diagonal 
in the basis of Eq. (7) 

(Sl2 ■ ^3) = ^ [S{S + 1) - Si2(Si2 + 1) - 53(53 + 1)] 

(fis ■ h) = \ [T{T + 1) - ti2(ti2 + 1) - ^3(^3 + 1)] • (13) 

To complete the evaluation, we need the matrix elements of the exponential. These can be 
obtained in analytic form 

j-QO 

InAo^ = e-" ^ [Rn,i{r)fr^dr . (14) 

Jo 

The results are straightforward. Here, by way of example, we quote the result for I = n—1 

In,l=n-l{a) = i—^r^' . (15) 
2r m 

Combining all pieces together, the Hamiltonian is 

2 

H = Eo + ^-- + Pr + iB + C6o)Ss,„i 
2m r 

+ (-l)'+i2Ae-"'-[sr2 ■ s-3 + ti2 ■ h + 2si2 ■ si fi2 ■ hi (16) 

where 5q stands in short for 5n,i^«,o- This Hamiltonian is characterized by the parameters 
r, /3, Eq, a, a, B, C, which are chosen by comparing with experimental data. Since all 
contributions are given in explicit analytic form, the determination of the parameters is 
straigthforward. The procedure that we use is the following: 

(i) the parameters r^m, ^ and Eq are determined from the location of the lowest state 
for each orbital angular momentum; we obtain r^m = 1546 MeV, — = 5 MeV and En = 
1706 MeV; 

(ii) the parameters B and C are determined by the splitting between the state S12 = 0, ti2 = 
and the average of the states S12 = 1, ti2 = 1- We find B = 300 MeV and C = 400 MeV. 
The value B = 300 MeV is consistent with earlier estimates, B = 250 MeV, arising from 
an evaluation of the color magnetic interaction ; 
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FIG. 2: Comparison between the calculated masses (black lines) of non-strange baryon resonances 
up to 2 GeV and the experimental masses (grey boxes) from PDG 22\ for the 4* and 3* resonances. 

(iii) the parameters A and a are determined from the spUtting of the multiplet within 
= 1, tu = 1- We find A = 205 MeV, ^ = 0.30. 

In Fig. 2 and in Table II the results of the model are compared with the experimental 
data 22]. It can be seen that the quark-diquark model with the specific interaction, Eq.(16), 
provides a good description of the masses of 4* and 3* resonances. The quality of this 
description is similar to that of the usual three-quark model in its various forms [2^, [2^, 
and it can be observed that the predicted value for the mass of the Roper resonance is 
higher than the experimental data. However, in the quark-diquark model, some degrees 
of freedom are frozen. There are therefore far fewer missing resonances, and in particular 
no missing resonances in the lower part of the spectrum. On the contrary, the problem of 
missing resonances plagues all models with three constituent quarks. 

It will also be noted that the clustering of states expected by the quark-diquark model is 
particularly evident for = 1~ and = 2"^, at both the qualitative and quantitative levels. 
This clustering appears both in the nucleon, N, and in the A. The clustering of states is 
another feature that is difficult to obtain in a three quark description. In the quark-diquark 
model, however, it is obtained automatically. 

For the model described here, in which all interactions in addition to the Coulomb-like 
interaction are treated in perturbation theory, elastic and transition form factors can all be 
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TABLE II: Comparison between the calculated masses of non-strange baryon resonances up to 2 
GeV, Mcaic, and the experimental masses 2j], M^xp-, for the 4* and 3* resonances. The masses 
are given in MeV. 
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calculated analytically. The scalar matrix elements are of the type 

j d\Rn'A^)yi',nA^)e~''^''Rn,i{^)yiA^) ■ (17) 
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These integrals, denoted by [/„ i,n'i'{k)Smm' are straightforward and Table III shows the 
corresponding results for transitions from the ground state with quantum numbers, n = 1, 
= 0+ to a state with n, /. It is observed that the elastic form factor is 



1 



[1 + W)2 ' 



(18) 



with a 



2rm ' 



In addition to having a power-law behavior with momentum transfer, k, 



typical of Coulomb-like interactions, this form factor has precisely the power dependence 
observed experimentally. Thus the quark-diquark model presented here has the further 
advantage of producing in first approximation an elastic form factor in agreement with 
experimental data. Fig. 3. All form factors depend on the scale a. To determine the scale 
a the r.m.s. radius can be calculated by using the ground state wave functions, < >= 



12a^. This calculated value is then fitted to the experimental value < > 



exp 



0.74(1) fm^ [22]. The resuhing value is a = 0.25 fm. Since the parameter r is determined 
from the spectrum, one obtains the reduced mass m = 102 MeV. This value is somewhat 
lower than the naive expectation m = 200 MeV, obtained by assuming the quark mass to 
be 300 MeV and the diquark mass to be 600 MeV. The results shown in Table III should 
be compared with the analogous results in the three quark model, as for example reported 
in Table IX of Ref. |2^ 
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FIG. 3: The electric form factor of the proton. The dotted line is the result of the model (Eq.(18)) 
using a = 0.25 fm, the dot-dashed line corresponds to the dipole fit, a = 0.23 fm. The experimental 
data are from Ref. 



I le Qc 

3- 



In order to calculate the magnetic elastic form factors and the helicity amplitudes, other 
matrix elements are also needed to be calculated; this program will be completed in another 
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TABLE III: The scalar form factors of Eq.(17) for transitions to final states labelled by the quantum 
numbers n, l^, where P is the parity. The initial state is n = 1,1^ = 0"*" and a = 
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article (not least because a relativistic version of the model is required for a good calculation 
of form factors), since here we are interested only to explore the qualitative features of the 
model and of its results. In particular the quark-diquark model presented here produces 
the phenomenon of stretching, which is at the basis of the Regge behavior of hadrons. 
The transition radii increase with n and as one can see from Table III, or by evaluating 
< i""^ >n,i=n~i= (2n + 2)(2n + l) ^jf^2 ■ In other words, hadrons swell as the angular 
momentum increases. 

In this article, we present a simple quark-diquark model with a specific direct plus ex- 
change interaction. This model reproduces the spectrum just as well as conventional three- 
quark models. However, it has far fewer missing resonances than the usual models. Most 
importantly, the model produces form factors with power-law behavior as a function of mo- 
mentum transfer, in agreement with experimental data. Finally, it shows the phenomenon 
of stretching which is at the basis of the Regge behaviour of hadrons. 

One may wonder whether there is a unique spectral signature for quark-diquark models. 
One of these signatures is the detection of 1"^ states which are antisymmetric in all three 
quarks. These states, originating from the omitted diquark representation 15 of SUsf{6) are 
not present in the quark-diquark model and occur (at different masses) in all models with 
three quarks. These missing states may, however, be very difficult to detect since they are 
decoupled and cannot be excited with electrons or photons. To excite these states, strongly 
interacting particles are needed, for example {p,p') with spin transfer. Another possibility 
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is to study whether or not the mixed symmetry states (which are doubly degenerate in the 

model) are in fact simply degenerate (as in the quark-diquark model). 

The work presented here can be expanded in several directions: (1) to include strange 
baryons; this expansion is straightforward and requires no further assumption; (2) to study 
multiquark states; (3) to include relativistic corrections. 

For transparency, here we have presented results in which all interactions, in addition to 
the Coulomb-like force, are treated in perturbation theory. Numerical diagonalization shows 
that with the parameters of this article, the perturbation treatment is valid (at least for the 
low-lying states). The complete numerical results will be the subject of a subsequent paper. 

An aspect that has not been discussed here is how to derive the quark-diquark model 
from microscopy. This aspect has been the subject of many investigations in the similar 
problems of condensed matter and nuclear physics. The usual argument in hadronic physics 
is that diquark correlations arise from the spin-spin interaction originating from one gluon 
exchange that lowers the scalar diquark relative to the vector diquark. Some interes ting 
results are now available and show that instanton interactions favor diquark clustering 14| . 
and that by using an Instanton Liquid Model these generate a deeply bound scalar antitriplet 
diquark not point-like In this respect instantons can provide a microscopic dynamical 
mechanism in terms of non perturbative QCD interactions. In this article, we have merely 
tried to describe many data by means of a phenomenological approach; the objective of a 
subsequent paper will be to understand what kind of microscopic QCD mechanisms we are 
trying to mimic in this oversimplified form. 
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